Decay of Quantum Correlations in Atom Optics Billiards with Chaotic and Mixed 

Dynamics 
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We perform echo spectroscopy on ultra cold atoms in atom optics billiards, to study their quantum 
dynamics. The detuning of the trapping laser is used to change the "perturbation", which causes 
a decay in the echo coherence. Two different regimes are observed: First, a perturbative regime in 
which the decay of echo coherence is non-monotonic and partial revivals of coherence are observed. 
These revivals are more pronounced in traps with mixed dynamics as compared to traps where 
the dynamics is fully chaotic. Next, for stronger perturbations, the decay becomes monotonic and 
independent of the strength of the perturbation. In this regime no clear distinction can be made 
between chaotic traps and traps with mixed dynamics. 



The decoherence of superposition states, as they inter- 
act with the environment |l| , is one of the main obstacles 
in realizing quantum information processing schemes . 
Decoherence or dephasing represent loss of information, 
quantified by a decrease in the "fidelity" , which typically 
measures the overlap between a desired state and the ac- 
tual output of a real system. 

In the field of classical chaos exponential sensitivity 
to initial conditions (non-vanishing Lyapunov exponent) 
plays a central role. However, the quantum analogue of 
evolving two initially slightly different wavefunctions in 
the same Hamiltonian yields that their overlap is con- 
stant in time, thereby not giving any information on the 
dynamics of the system. The fingerprints of the classi- 
cally chaotic dynamics in the quantum regime (Quan- 
tum chaos) is a rich field of study 3] . One approach is to 
study the decay of the overlap between two initially iden- 
tical wavefunctions evolved in slightly different Hamilto- 
nians. The fidelity then denotes the overlap between a 
state evolved by a Hamiltonian Hf with the same state 
evolved by a slightly perturbed Hamiltonian H± jj| . The 
fidelity is in this context often denoted the Loschmidt 
Echo, since it is equivalent to the overlap between an 
initial state, and the same state evolved forward in time 
in iff and then backwards in time in H± [||. The de- 
cay of fidelity in chaotic systems and its dependence on 
several parameters, has been the topic of intense theo- 
retical studies in recent years (see Ref. 0, 0, @ and ref- 
erences therein). Nevertheless, experimental studies of 
chaotic systems are still lacking (mostly due to the dif- 
ficulty of preparing highly-excited pure quantum states) 
and so are both theoretical and experimental investiga- 
tions of systems with mixed dynamics. Of special inter- 
est in quantum dynamical studies are the "cross-over" 
regimes through which the system goes by changing its 
parameters from those values for which a quantum de- 
scription is required, to those that allow classical models 
to be used. Understanding these regimes can shed light 
on the relation between classical and quantum properties 



of the system. 

Ultra cold atoms have been used in the past to ex- 
perimentally study both quantum and classical dynam- 
ics. Quantum dynamics have been studied in driven ID 
systems, where a broad variety of phenomena such as 
dynamical localization Q, dynamical tunneling and 
quantum accelerator modes Q have been demonstrated. 
Classical dynamics has been studied in atom-optics bil- 
liards [id. in which regular, chaotic and mixed dy- 
namics were observed. 

In this letter we use microwave (MW) "echo spec- 
troscopy" 01 to experimentally measure the dephas- 
ing and quantum dynamics of ultra cold 85 Rb atoms 
trapped in atom-optics billiards, with underlying chaotic 
or mixed classical dynamics. Echo spectroscopy mea- 
sures the overlap between two initially identical states 
evolved in slightly different Hamiltonians (and is there- 
fore closely related to the fidelity |l3j), but it overcomes 
the need for pure quantum states and allows the use of 
thermal ensembles for the study of quantum dynamics 
[l2|. We demonstrate that the decay of the echo co- 
herence displays qualitative different behavior for differ- 
ent perturbation strengths. Two distinct regimes are ob- 
served. First, a "perturbative" regime where the decay 
depends on the perturbation strength and where partial 
coherence (wavepacket) revivals are observed even when 
the underlying classical dynamics is chaotic. This indi- 
cates that the partial revivals observed for a nearly har- 
monic trap |12| are in fact a generic feature of trapped 
atoms, as predicted in [l3|. For stronger perturbations 
we observe a crossover to a regime where the decay of 
the echo coherence is independent of the perturbation 
strength. In this regime the decay is monotonic and no 
revivals are observed. For traps where the classical mo- 
tion exhibits a mixed phase space (i.e. stable "islands" in 
a chaotic " sea" ) we observe more pronounced revivals in 
the perturbative regime, due to the periodic classical mo- 
tion in the islands. However, a perturbation-independent 
regime exists also for these systems, and then the decay is 
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essentially indistinguishable from that of traps with fully 
chaotic dynamics. 

We use 85 Rb atoms in a coherent superposition of 
their two magnetic-insensitive hypcrfinc Zccmann ground 
states. These two levels, \5Si/ 2 , F = 2,niF = 0) de- 
noted ||), and 5Si/ 2 , F — 3, m F — 0} denoted |f), are 
separated by the energy splitting Ehf = ^hf with 
luhf — 2irx 3.036 s . The atoms are trapped in a dipole 
potential formed by a linearly polarized laser close to the 
5Si/2 — > 5P 3 / 2 transition. The dipole potential is in- 
versely proportional to the trap laser detuning , hence 
it is slightly different for atoms in ||) and |). The ex- 
ternal (center of mass) potential depends on the internal 
(spin) state, hence the internal and external degrees of 
freedom can not be separated, and the entire Hamiltonian 
(neglecting interactions between the atoms) is written as: 

h = ffiU)ai + (H T + £Hf)|T>(T| 
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where V± and Vj are the external potentials for atoms 
in states ||) and ||), respectively. These potentials in- 
clude the gravitational potential, equal for both states, 
and the dipole potential, which can be written as V^.j. and 
Vd,-\ = (1 + e)Vd,i, wnere e = ^hf/Al is the "perturba- 
tion strength" typically 10~ 3 — in our experiments 

The eigenenergies of this Hamiltonian consists of two 
manifolds (belonging to ||) and ||)) separated in energy 
by Ehf ■ The atoms are initially prepared in their inter- 
nal ground state ||) and their total wavefunction can be 
written as '5 = ||) ® ip, where ip represents the center of 
mass part of their wavefunction. The echo sequence con- 
sists of three short and strong MW pulses, each of which 
changes the internal state of the atoms, while leaving 
the center of mass part of the wavefunction unchanged 
[l2|. First a 7r/2-pulse puts the atoms into a coherent 
superposition of ||) and ||). After a time t a 7r-pulse in- 
verts the populations and after another time r the atoms 
are irradiated by a second 7r/2-pulse. The populations 
of ||) and |t) are then measured. If we start with an 
eigenstate |nj.) of Hi then the population of ||) after 
the echo pulse sequence is 0: Pj = ^ [1 — Re (F echo )], 
where F echo = {n i \e lH i T e iH ^ T e~ iH i T e~ iH ^ T \n i ) is de- 
noted the "echo amplitude" (we omit H to simplify the 
equations). F ec ho — 1 indicates perfect coherence and 
yields Pf = and F ec ho — yielding Pj = 0.5 indi- 
cates complete loss of coherence. If e = then internal 
and center of mass motion degrees of freedom decouple 
and Fecho = 1 for all times. When considering eigen- 
states the echo amplitude can be written as a time cor- 
relation function P ec / lo = e ZUJ "' lT (ip n (t = 0) ip n (t = r)) 
where \<p n (t = 0)) = e~ iH i T |n x ) and \ip n (t = r)) = 
e - lH i T \ l p n [t = 0)). Therefore, the decay of the echo co- 
herence corresponds to a decay of quantum correlations 
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FIG. 1: (Color online) CCD images of edge billiards with 
a = 52.5° (A) and a = 31° (C), used in our experiments and 
their corresponding calculated Poincare surface of sections, 
respectively, indicating chaotic dynamics (B) and mixed dy- 
namics with large islands of stability (D). 



due to dynamics in the trap. 

In our experiment 85 Rb atoms are loaded into a far 
off resonance optical trap, cooled to a temperature of 20 
fj,K, and optically pumped into the F=2 hyperfine state. 
By changing the detuning of the trap laser, and simulta- 
neously adjusting its power, the perturbation strength is 
controlled. After the MW echo pulses Pj, the population 
of state |t), is measured using fluorescence detection, and 
the signal is normalized to Pj- after a short 7r-pulse. 

The trap is a light-sheet wedge billiard, made from 
two crossed blue detuned light sheets defining the bil- 
liard walls and where gravity confines the atoms in the 
vertical direction [l4|. The light sheets have (1/e 2 ) di- 
mensions of 20 x 250 /im, and by the use of cylindrical 
lenses mounted on rotational stages, the wedge angle are 
adjusted in order to control the classical dynamics. The 
very elongated shape of the trap allows us to consider 
only the transverse motion and neglect the longitudinal 
one, which has a timescale much longer than the experi- 
ment time. The temperature of the atoms is much larger 
than the mean level spacing in the trap, and the atoms 
typically occupy many (up to ~ 10 8 ) states in the trap. 
The measured echo signal is the ensemble average of all 
of them. 

The structure of phase space in the wedge billiard can 
be tuned from stability to chaos by varying the vertex 
half-angle, a 0. For a < 45° phase space is mixed, and 
the size of the stable islands oscillates as a function of a, 
a behavior dubbed "breathing chaos" . For a > 45° the 
system is fully chaotic. Figures IB and ID present the 
Poincare surface of section for the two wedge billiards of 
Fig. 1A and 1C, respectively, used in our experiments, 
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calculated by numerical integration of the classical tra- 
jectories in the measured billiard potentials. As seen, 
the classical dynamics is indeed almost fully chaotic and 
mixed for a = 52.5° and a = 31°, respectively, as pre- 
dicted and previo usly measured [Io| . 

As shown in Ref. 13], the echo amplitude for small 
perturbations can be written as: 

Fech (K),t) ^4^] exp(-iw m r) |(m T | n x )| 2 - 

^2 exp(-2iw m r) |(m T | n x }| 2 + |(n T | n x )| 6 , (2) 

This means that the perturbative echo signal for a mo- 
noenergetic ensemble of atoms can be expressed as a func- 
tion of the local density of states (LDOS). The LDOS 
denotes the local average of the absolute value squared 
of the matrix elements of the transformation matrix from 
eigenstates of Hi to eigenstates of Pff . Formally, it is sim- 
ply \(mi \ni) | as a function of — averaged over a 
ensemble of neighboring with approximately the same 
energy. The width over which the LDOS is nonvanishing 
is denoted the "bandwidth" , and if it is large, we expect 
a rapid decay of the echo coherence. In the perturbative 
regime the LDOS displays system specific features, de- 
spite the fact that the underlying classical dynamics is 
chaotic. In particular it is evident from the semiclassi- 
cal calculations of that for atom optics billiards, 

where the inherent perturbation is localized on the bil- 
liard walls, the LDOS will have pronounced peaks for 
77i | — 71 J corresponding to E n — E m = h/ru, where tu is 
the typical time between encounters with the wall (~ 15 
ms in our experiments). This means that the echo am- 
plitude will show partial revivals for r = Tu fl3| . 

In Fig. |2 the decay of the echo signal for different per- 
turbations is presented for a wedge with a = 52.5° (see 
Fig. where the dynamics is almost fully chaotic, as 
shown in Fig. ^3. For small perturbations a nonmono- 
tonic decay is seen with a partial revival of correlations 
for t ~ tu as predicted in [Lj . The revivals are seen de- 
spite the fact that due to the "high" temperature a large 
band of energies is occupied in the trap. Partial revivals 
of correlations in traps where the dynamics is separable 
are generally expected at time scales equivalent to the 
ID level spacing, but their surprising observation here, 
for a trap with chaotic motion, demonstrates that they 
are a far more widespread phenomena. 

For larger perturbations a crossover to a regime where 
the decay is monotonic is observed. In this regime the 
decay is independent of perturbation strength. This is 
evident from Fig. in which Pj (t = 0.0025s) is plot- 
ted as a function of perturbation strength. It is seen 
that Pf initially grows with perturbation strength but 
for e > 0.004, Pj is roughly constant. The perturbation 
independent regime of the echo decay is associated with 
the nonuniversal regime of (ltij . where the LDOS was 
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FIG. 2: Echo signal for a light sheet wedge with chaotic 
classical dynamics, for different perturbation strengths (trap 
laser wavelength tuned from A = 775.9 to A = 779.7). ♦: 
e = 1.44 x 10" 3 , V: e = 1.90 x 10" 3 , A: t = 2.43 x 10~ 3 , o: 
e = 3.80 x 10~ 3 and ■: e = 1.52 x 10~ 2 . For small perturba- 
tions a nonmonotonic decay with revivals around r = tu is 
seen, whereas for large ones a monotonic decay is observed. 




FIG. 3: Echo signal for a time between pulses of 2.5 ms, 
as a function of perturbation strength. For a small perturba- 
tion the decay depends strongly on the perturbation strength, 
whereas for e > 4 x 10 -3 the decay is almost independent of 
perturbation strength. Dashed line: Pj (at 2.5 ms) predicted 
by the simple classical model described in the text. Inset: 
Solid line: Pf calculated by classical model. ■: Pf measured 
for a perturbation of e = 1.52 x 10~ 2 . 



found to be perturbation independent. In this regime the 
perturbation is large enough so the overlap of equivalent 
eigenstates is small, and this indicates that the effects of 
quantization of the trap levels should not play a role and 
a classical description might be possible. Since the echo 
amplitude can be viewed as a propagator, and in our sys- 
tem the thermal de-Broglie wavelength is much smaller 
than the billiard's dimensions, it is possible to use a semi- 
classical propagator [3( . It is then seen that the classical 
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FIG. 4: Echo signal for a light sheet wedge where the classical 
phase-space is mixed. •: e = 1.44 x 10 -3 . □: e = 1.52 x 10~ 2 . 
For the small perturbation a revival that is more pronounced 
than the one in the trap with classical chaotic dynamics is 
observed, whereas for large perturbations the decay can not be 
distinguished from the decay in a trap with chaotic dynamics. 



trajectories contributing to the ensemble average of the 
echo amplitude, are those that after evolving forward in 
time in H-\ and H±, and then backwards in time in 
and Hi, return to the vicinity of their initial position. 
These trajectories we divide into two types: those that 
during this propagation have hit the wall, and those that 
did not. Since Hy and H± are highly different mainly in 
the vicinity of the wall, then the action integral along the 
first type of trajectories will yield a very large phase and 
the contribution from these trajectories to the ensemble 
average of the echo amplitude will average out. However 
since the second type of trajectories does not feel the dif- 
ference between the potentials it will retrace its forward 
propagation backwards in time causing the action inte- 
gral to vanish. These trajectories will give a perfect con- 
tribution. Therefore the echo amplitude in this regime 
simply measures the probability that the particles have 
not yet hit the wall, and this is a classical quantity. The 
dashed line in Fig. [3] describes the classical estimation of 
Pf (t = 0.0025s) calculated assuming an idealized hard 
wall wedge populated with a thermal ensemble of atoms 
at a temperature of 20 fj,K, clipped at an energy equal to 
the depth of the trap. In the inset of Fig. [31 the classical 
calculation is shown together with the measured echo de- 
cay for a perturbation of e = 1.5 x 10~ 2 , and reasonable 
agreement is seen despite the extreme simplicity of the 
above model. 

Next we consider billiards with mixed dynamics. The 
semiclassical perturbative calculations of 0, 0] indi- 
cate that classical periodic motion as in elliptic islands 
will lead to narrow peaks in the LDOS, thereby yield- 
ing more pronounced revivals compared with the chaotic 
motion discussed above. This is observed in Fig. 0] 
where the decay of the echo signal is seen for the mixed- 



dynamics billiard with a = 31°, shown in Fig. ^ The re- 
vival for small perturbation is deeper than for the chaotic 
motion of Fig. and the reminiscence of a second re- 
vival for t ~ 2tu can be seen indicating that the peak in 
the LDOS is narrower. Since the revivals are more pro- 
nounced in the case of mixed phase space they "survive" 
to slightly larger perturbation strength. However, in the 
perturbation-independent regime also shown in Fig. 
the decay is essentially indistinguishable from that of the 
chaotic billiard, in agreement with the classical model 
given above. 

In summary, we studied the decay of quantum cor- 
relations in atom optics billiards in which the classical 
dynamics is chaotic or mixed. We observed two dis- 
tinct regimes for the perturbation strength in which the 
decay was qualitatively different. In the perturbative 
regime the decay was non-monotonic, with revivals at 
times corresponding to the typical time between bounces 
from the wall. The revivals were more pronounced in 
traps with mixed phase-space as compared with traps 
where the dynamics is almost fully chaotic. In the per- 
turbative regime the decay rate increased with pertur- 
bation strength. However, in the perturbation indepen- 
dent regime the decay was monotonic and independent of 
perturbation strength. No clear distinction between the 
decay in traps with mixed and chaotic dynamics could 
be made in the perturbation independent regime. In this 
regime the echo amplitude measures the classical proba- 
bility that an atom have not yet hit the wall. 
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